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The occurrence of antiferromagnetism in k-(ET)2X can 
be understood within an effective i-filled band with 
dimers of ET molecules containing one hole each. We 
argue that while this effective model can describe the 
presence of antiferromagnetism, a complete description 
for these materials requires the correct carrier density of 
one-half per molecule. For dimerized and strongly frus- 
trated i -filled lattices we show that a singlet-paired state 



coexisting with charge ordering occurs that we have 
termed the Paired Electron Crystal (PEC). Here we in- 
vestigate the i -filled model on a dimerized lattice, show- 
ing regions where AFM, PEC, and the Wigner-crystal 
occur. We point out the need to go beyond quantum spin 
liquid concepts for highly frustrated materials such as k- 
(ET)2Cu2(CN)3 and /3'-EtMe3Sb[Pd(dmit)2]2 which we 
believe are PECs at low temperatures. 
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1 Introduction A challenging problem in condensed 
matter physics is to describe systems where both electron- 
electron (e-e) interactions and lattice frustration are signif- 
icant. The possibility that on a frustrated lattice a quantum 
spin system may remain disordered even at low temper- 
ature forming a quantum spin liquid (QSL), has spurred 
much theoretical and experimental work [1]. Two candi- 
date QSL materials, k-(ET)2Cu2(CN)3 (hereafter k-CN) 
and ^'-EtMe3Sb[Pd(dmit)2]2 (hereafter dmit-Sb), have 
been identified from within the group of charge transfer 
solid (CTS) organic superconductors [2]. In addition to 
these particular two compounds displaying QSL-Iike be- 
havior, the K-(ET)2 and /?'-X[Pd(dmit)2]2 series are of 
wide interest because they display other exotic phases 
including antiferromagnetism (AFM), a so-called valence- 
bond solid (VBS) phase, and superconductivity (SC) [3, 4]. 

Structurally both k-CN and dmit-Sb consist of two- 
dimensional layers of dimers arranged in an anisotropic 
triangular lattice. With an average of 0.5 carriers per 
molecule this has led to the common description of both 



within an effective i-filled band where each dimer is re- 
placed by a single site. In the presence of strong on-site e-e 
interactions one obtains a quantum spin model where each 
spin represents a pair of molecules [2]. 

The i -filled Hubbard model on an anisotropic triangu- 
lar lattice is frequently taken as a minimal electronic model 
for these materials [5]. This model contains a single frus- 
trating bond t' per square plaquette. In the case of weak 
anisotropy (small t' /t), the ground state has Neel antifer- 
romagnetic order As t' /t is increased, frustration of the 
AFM leads to a paramagnetic metallic (PM) state. Based 
on mean-field calculations many investigators have sug- 
gested that a region of (1^2 _y2 SC mediated by spin fluctua- 
tions emerges near the boundary between the AFM and PM 
phases. This proposal ignores other CTS where the super- 
conducting state is proximate not to the AFM but to charge 
ordered or VBS phases. More importantly, calculations be- 
yond the mean-field approximation do not find SC within 
the i -filled model for any U and t'/t [6,7]. Our earlier 
exact calculations were for a finite 4x4 cluster [6]. More 
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Figure 1 (a) Dimer lattice considered here. Double lines 
correspond to rigid bonds with hopping matrix elements 
t + St. Dashed and double dashed lines have hopping ma- 
trix elements of t — St in the absence of lattice distortion. 
Electron-phonon coupling is included on inter-dimer bonds 
in X and y directions, (b) PEC state. Here filled (empty) 
sites have charge density of 0.5+5 (0.5-5). (c) The WC-SG 
state. In (b) and (c), heavy lines indicate singlet bonds. 




Figure 2 Charge disproportionation An in the dimer lat- 
tice versus t'. Parameters are U ~ 6, Vx = Vy = 1, 
V' = 0, a = 1.0, 13 = 0.1, Ka = Kp = 2. Diamonds, 
circles, and squares are for x-axis dimerization of St =0.2, 
0.3, and 0.4, respectively. Lines are guides to the eye. 



recently, we have adopted the Path Integral Renormaliza- 
tion Group method [8] to perform calculations for much 
larger lattices confirming the absence of SC [9]. 

The absence of SC within the i -filled effective model 
for the CTS suggests that one must return to the underlying 
I -filled lattice in order to fully understand the electronic 
properties. We have recently shown that in the ^-filled 
band, strong e-e correlations and lattice frustration lead to 
a transition from AFM to a paired charge-ordered state, 
the Paired Electron Crystal (PEC) [10, 11]. The PEC can 
be visualized as the ^-filled equivalent of a VBS state in 
a quantum spin model. In the i-filled picture, each single 
spin of the spin model must be replaced by a pair (dimer) 
of sites. Once singlet bonds are formed between dimers, 
the charge occupation within each dimer necessarily be- 
comes nonuniform — unlike the ^ -filled band, at i -filling 
charge and spin degrees of freedom are linked. Interest- 
ingly, recent experiments (discussed below in Section 3) 
have found evidence for charge disproportionation in the 
optical spectra and the frequency dependence of the dielec- 
tric response at low temperature. 

In references [10,11] we considered a |; -filled an- 
isotropic triangular lattice with Coulomb and electron- 
phonon (e-p) interactions where along any one direction 
(x, y or X + y) all bonds were taken to be uniform in the 
absence of instabilities. In order to simulate the dimerized 
K or dmit lattice, we have subsequently repeated these 
calculations for a lattice in which the sites have formed 
rigid dimer units along one direction. We report these re- 
sults here. The organization of the paper is as follows. In 
Section 2 we describe calculations showing that the PEC 
state occurs also within the lattice of coupled dimers. In 
Section 3 we discuss recent experiments for the proposed 
CTS QSL materials, that suggest that at the lowest tem- 
peratures a I; -filled model is essential to correctly describe 
their electronic properties. 



2 Dimer lattice calculations We consider the fol- 
lowing Hamiltonian, 

H = - J2 ^-(1 + "-^y )^^.- + ^ E ^"^y 

i i 

In Eq. 1, Bij = J^ai'^l a'^j,'^ + H.c.) is the kinetic en- 
ergy operator between sites i and j. The index i' refers 
to the lattice directions x, y, and x + y (see Fig. 1). a^, 
is the inter-site e-p coupling for direction with corre- 
sponding spring constant A'^. /3 is the onsite (Holstein) 
e-p coupling with corresponding spring constant Kp. U 
and V are respectively the onsite and nearest-neighbor 
Coulomb interactions. We have previously considered the 
case where the intersite e-p coupHng applies to all t^ and ty 
bonds [10, II]. Here we will consider a lattice consisting 
of rigidly bound dimers connected by deformable bonds, 
with frustrating bonds t' = t^+y across each plaquette of 
the lattice as shown in Fig. 1(a). In this lattice t^ varies 
as tx + (— l)'^(5i, where ix is a site index in the x direc- 
tion. The bond modulation induced by the e-p coupling 
^iT.iT+i is set to zero for the strong intra-dimer bonds, 
and is allowed to vary self-consistently for the inter-dimer 
bonds along the x-direction. E-p coupling along the y axis 
is not restricted. As in reference [11], we take ax = ay, 
Kx = Ky = Kp = 2, tx = ty = t, and measure energies 
in units of t. We solve Eq. 1 exactly on a 4 x 4 lattice. 

As in the lattice without rigid dimerization [10,11], 
when the strength of the frustrating bond t' is increased we 
find a transition from AFM (Fig. 1(a) with uniform charge 
density on each site) to a charge-disproportionated PEC 
state (Fig. 1(b)). Fig. 2 shows the charge disproportionation 
An, defined here as the charge density difference between 
sites with large and small charge densities, as a function of 
t'. For t' > t'^. An becomes non-zero indicating charge 
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Figure 3 Variation of (a) charge disproportionation. An, 
and (b) bond order between two successive bonds in the 
X + y direction with V. Nearest-neighbor Coulomb inter- 
actions of the form = Vy = V and V' = Q are assumed. 
Other parameters are U = 6, t' = 0.6, a — 1.0, /? = 0.1, 
Ka = Kp = 2. Lines are guides to the eye. 



disproportionation within each dimer. Simultaneous loss 
of antiferromagnetic order and the opening of a spin gap 
due to nearest-neighbor singlet formation as indicated in 
Fig. 1(b) is seen in the bond orders and spin-spin correla- 
tions [10, 11] (not shown here). Assuming equal e-p cou- 
pling strength, we find that stronger dimerization (larger 
6t) gives the PEC state at a weaker level of lattice frustra- 
tion. However, as is seen in Fig. 2, for strong dimerization. 
An remains relatively small until t' is significantly larger 
than t'^. This may have experimental implications. 

When nearest-neighbor Coulomb interactions Vij in 
Eq. 1 are stronger, the checkerboard-pattern Wigner Crys- 
tal (WC) state shown in Fig. 1(c) can become more stable 
than the PEC of Fig. 1(b). Here we choose Vij to be of 
the form ^ Vy = V but V = 0. Other choices of 
Vij can change the stability of the WC significantly — for 
example taking V' > suppresses the checkerboard WC 
shown in Fig. 1(c) and stabilizes the PEC. [11]. As shown 
in Fig. 3(a), Ari in the WC state is considerably larger than 
in the PEC state. The charge and bond distortion of the WC 
state is shown in Fig. 1(c). Because of the rigid dimers in 
this lattice, the bond strength alternates strong-weak-strong 
along the x + y (<') direction. The bond alternation in the 
one dimensional chains along the x + y direction gives this 
phase a spin gap, denoted as "WC-SG" in reference [11]. 
In Fig. 4 we show the t' — V ground state phase diagram. 

Summarizing our numerical results for the case the lat- 
tice of rigid dimers; (i) As in the lattice with unconstrained 
tx bonds we find an AFM-PEC transition at a critical frus- 
tration t[. Stronger dimerization leads to PEC formation 
over a larger region of phase space, however with a weaker 
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Figure 4 Ground state phase diagram of the dimer model 
for 6t = 0.3, U^6,V^=Vy^V, and V = 0. Lines are 
guides to the eye. 



An close to t[. (ii) When nearest-neighbor Coulomb inter- 
actions are significant a WC phase also occurs. In the case 
of a rigid dimer lattice this WC phase has a much greater 
tendency towards spin gap. 

3 CIS spin liquid candidates In this section we dis- 
cuss K-CN and dmit-Sb in detail. SC occurs in both: in k- 
CN underpressure [12], and in other dmit materials [13]. 
For both materials a number of recent experiments indicate 
transitions to broken-symmetry states at low temperature. 

3.1 /«-(ET)2CU2(CN)3 In K-CN the estimated Heisen- 
berg exchange integral between nearest-neighbor dimers 
is J ~ 220-250 K [14]. Despite the significant magnetic 
interaction between dimers ^H NMR experiments find no 
long range magnetic order to mK temperatures [14]. How- 
ever, a number of experiments suggest that a second order 
phase transition occurs at '-^ 6 K [15-20]. Below 6K the 
^•^C NMR relaxation rate 1/Ti decreases, suggesting a spin 
gapped ground state [15]. A decrease in the spin fluctu- 
ation rate is also seen below 6K in zero-field /iSR [20]. 
Field-dependent /iSR finds a critical field He above which 
the linewidth increases rapidly and estimate the spin gap as 
3.5 mK [20]. Precision measurement of lattice expansion 
finds strong lattice effects suggesting a second-order phase 
transition at 6K [19]. The dielectric response is frequency 
dependent below 60 K suggesting an antiferroelectric or- 
dering of dipoles below 6 K, requiring unequal charge den- 
sities within each dimer [18]. One of the most puzzling 
aspects of is that while certain measurements (^'^C NMR, 
yuSR, and thermal conductivity [17]) indicate a gapped 
state, the specific heat is linear in T at low temperature con- 
sistent with a gapless energy spectrum [16]. 

3.2 EtMe3Sb[Pd(dmit)2]2 dmit-Sb is similar in 
many respects to k-CN: in terms of dimer units, the frustra- 
tion t' /t and effective J for the two materials are similar. 
No apparent magnetic ordering occurs down to mK tem- 
perature [21]. Frequency-dependent dielectric response 
similar to that of k-CN is found at low temperatures [22]. 
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In the related material EtMe3P[Pd(dmit)2]2 described as 
having a VBS ground state [23], peak splitting below 25 
K in the infrared optical response indicates that the VBS 
state is charge ordered [24] . The charge order patterns in 
EtMesP are in complete agreement with those in the PEC 
[10, 11]. In dmit-Sb, the same IR peaks are not split but 
broadened, prompting the investigators to suggest fluc- 
tuating charge and bond order [24]. Similar to the 6 K 
transition in k-CN, a transition to a broken-symmetry state 
appears to occur in dmit-Sb at T '-1 K. The ^^C NMR 
relaxation rate 1/Ti decreases below 1 K, but with a dif- 
ferent temperature dependence than in k-CN [15,21,25]. 
Thermal conductivity measurements in weak magnetic 
fields suggest gapless mobile excitations, but the steep in- 
crease of the thermal conductivity when H >2 T shows 
the presence of gapped spin excitations [26]. 

3.3 Conclusion Within the effective i-filled model 
little difference would be expected between k-CN and 
dmit-Sb as both have similar values of inter-dimer t' /t. 
While both materials display a mix of gapped and gap- 
less behavior, there are significant differences between 
the two such as the temperature dependence of the ^'^C 
l/Ti, which is oc ri in k-CN [15] but cx in dmit-Sb 
[21, 25]. These and other differences suggest that a i-filled 
description is essential at low temperature. We have sug- 
gested that the ground states of k-CN and dmit-Sb are 
PECs [10, 11]. In such a picture the state below 6 K in 
K-CN and 1 K in dmit-Sb would be a charge dispropor- 
tionated spin gap state with fluctuating short-range order 
/xSR gives clear evidence for a small but finite spin gap 
from the ground state in k-CN[20], in agreement with 
our prediction. A second transition at higher field[20] is, 
however, currently not understood. An alternate plausible 
model is the non-magnetic insulator phase found in the 
i-filled model, which however has no spin gap [27]. 

A paired state like the PEC can also explain the ob- 
served mix of gapped and gapless behavior. A near de- 
generacy of spin-singlet configurations can give a linear 
specific heat. This degeneracy is particularly expected in 
K-CN, where the lattice structure suggests many nearly 
equivalent ways to construct the spin-singlet bonds of the 
PEC [10]. On the other hand, magnetic excitations require 
breaking singlet bonds and will be gapped. A similar com- 
bination of gapless specific heat but gapped magnetic sus- 
ceptibility and conductivity is found in vanadium bronzes 
[28]. It was suggested singlet inter-site bipolarons could 
explain these features — tunneling of the carriers compos- 
ing the bipolarons between nearly-equivalent lattice sites 
gives a finite degeneracy and linear specific heat, while 
magnetic susceptibility is gapped because of the pairing 
gap of the bipolarons [28]. 

The suggestion that a paired state can describe the low 
temperature properties of k-CN and dmit-Sb also gives a 
possible mechanism of superconducting pairing [10,11, 
29,30]. In the picture, spin-singlet pairs formed in real 
space would form a superconducting paired-electron liquid 



[31]. Within a simplified model where pairs of sites in the 
i -filled model are replaced by single sites in an effective 
i-filled model with attractive {—U) interactions, a charge- 
order-to-SC transition occurs in the presence of lattice frus- 
tration [29]. We emphasize that our proposed mechanism 
is different from the mean-field theory of charge fluctua- 
tion driven SC [32]. Strong correlation, i-filling, and lat- 
tice frustration are also characteristic of unconventional SC 
in other families of materials, including the layered cobal- 
tates, spinels, and fullerides [30]. 
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